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e lom: ot *positiyvely cnharged particles that аге 
Generated on or near satellite surfaces have been made оп 
several spacecraft. This thesis postulates sputtering of 
the satellite surface due to ambient ion impact as the 
generating mechanism. Calculations are made using the 
Sigmund-Thompson sputtering theory to determine the response 
at the particle detectors. These calculations indicate that 
surface sputtering creates a sufficient flux to account for 
the observed phenomena. The NASA Charging Analyzer Program 
was run to determine the trajectories for actually observed 
particles. The calculated trajectories were determined to 
lead to the spacecraft surface, again indicating that 
surface emission was the source. Tie=Spucrteying u flüx as 
calculated was insufficient to cause any significant short- 


term damage to the spacecraft, beyond thin coating erosion. 
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1, INTRODUCTION 


A. OVERVIEW 

Spaeecraft charging is the development of a potential 
difference between a space vehicle and its plasma 
Коле. || lo OL slOnifticance for scientific purposes, 
since the potential may affect measurement of the ambient 
оре јез Of tie Environment, and for practical purposes, 
since it may lead to anomalous command signals and physical 
damage to the affected spacecraft. ШИ houen реса толо о 
Spacecraft charging had been made previously, the first 
МЕЧЕ Со of high level spacecraft charging were those for 
moelications Technology Satellite (ATS) 5 {Ref. 1}. Similar 
events were subsequently noted on ATS-6 and P78-2 (SCATHA) 
by various observers {Ref.2}. 

Generally, negative charging events are identified ру ап 
intense peak in the measured ion spectra at a specific 
energy, with an absence of ions of lesser energies. TRIS 
шшеп-=е peak 15 most clearly seen on the plot of the 
w ribution function versus energy, as will be displayed 
Mer in this work. This peak is known as the charging 
эг Such observations result from the acceleration of the 
ambient plasma ions into the detector by the potential drop 
between the plasma and the spacecraft. Occasional y; So 


fluxes are seen at energies below the peak energy, which 


would appear to violate conservation of energy вини 
collisionless “plasma” These ion distributions have been 
reported previously, and tentatively ascribed to sput ermi 


from spacecraft surfaces {Ref. 3}. They have been termed 


ry 


"spacecraft generated ons Бер > 

Analysis of this particular phenomenon is of рїе та 
interest, as these ions may contaminate the charging peak, 
leading to an underestimation of the satellite россии 
The contamination of the charging peak may be рак ош ан. 
important active experiments in low earth огр ,/50т 0 8 
ЭӘРАСЕБАВ- Additionally, since we believe the source to be 
sputtering, the flux of these particles is an indica p 
the damage to the satellite surface caused by the We yen 
environment. Also, the life expectancy of precision 
satellite surfaces, e.g. optical surfaces, and surface 
coatings may be directly affected by the sputtering tate 


appropriate to their environment. 


B. SCATHA SPACECRAFT 

Most of the experimental data and calculations performed 
in this thesis are based on the P78-2 (SCATHA) spacecrait. 
The construction, environment, and sensors ОР Спіс УӘ ЕЕ 
Will be discussed in the nekt several sections. 

1. o SpacecrareEMeoöne emre e Men 

The U.S. Air Force P78-2 spacecraft vas launchec mm" 

January 1979 to collect data for a study о: Space rum 


Claram AT High Altitudes юОСАТНД 111 а jora: 


NASA/Department of Defense Progran. The satellite body was 


а спаг1сај]ј in shape with a length and diameter of 
approximately 1.75 meters. Seven booms were deployed in 
шої С from the spacecraft to provide isolation Fox 


experiments from charging effects on the satellite surface. 
Figure 1 depicts the SCATHA spacecraft. The top of the 
Satellite in the figure is the forward end, and the bottom 
ако. The experiments and communications equipment for 
the satellite are primarily located on equipment decks in 
the center of the cylinder. The University of California, 
San Diego, charged particle detector was on the forward end, 
and the ion gun was on the aft end; both are described 
below. The outer cy mnd resurface was divided into three 
general areas; two solar arrays, one forward and one aft, 
and a bellyband between them to provide access panels. 
Additionally, sections of the surface were covered in 
various Materials, such as kevlar and mylar, for use in the 
experiments. {Ref. 4} 

The spacecraft was constructen GL special 
lightweight materials, as its orbit was at the maximum 
capability of the launch vehicle. mercat ral Lube of tie 
spacecraft was magnesium and the ecuipment decks were 
И лит. mem OlT arrays were aluminum core with a 
a glass (S102) outer face. The bellyband panels were 


covered with thermal paint and second surface mirrors, апа 
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SCATHA Satellite (Ref. 


Figure 1: 


acted as waste heat radiators. The forward cylinder end was 
coated with gold, except for some sample patches. {Ref. 4} 
2. Spacecraft Orbit and Environment 

The vehicle was inserted into a near-geosynchronous, 
eliptical orbit with a perigee of 5.5 Ra and an apogee of 
2 Ra. The orbit period маз 2325 hours, and the 
inclination to the equatorial plane was 7.8 degrees. The 
satellite axis was perpendicular to the earth-sun line, and 
the vehicle had a spin period of 59 seconds. The satellite 


encountered 40 day eclipse periods in both the fall and the 


spring. During these times, the spacecraft was shadowed for 
approximately 1 hour per day. Figure 2 gives a schematic 
веш спе Orbit. The view is from above the earth's north 


ole, $Ref.5] 

The environment of the SCATHA spacecraft consisted 
СЕ two main regimes, the plasmasphere and the plasmasheet. 
The nature of these regimes, particularly the ion 
composition, is important in determining the magnitude of 
the vehicle's interaction with the environment, for example, 
the sputtering yield. The boundary between these regimes 
(the plasmapause) is well defined, but varies with solar 
activity and solar wind parameters. The spatial regions are 
depicted in Figure 3. The view depicts the northern 
hemisphere. {Ref. 6} 

The innermost region is the plasmasphere. The 


density is generally greater than 105 т 7. and the 
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SCATHA Orbit {Ref.4} 


Figure 2.: 
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Regions of Space near Earth (Ref. 


Figure 3. 


temperature is less than 1 eV. The dominant ion 15 ПТ”) 
о“ and Не“ contributing approximately 20% ot the 1101 
density. The source of these ions is currently conSsi и 
to be the lono phere. (ReiS. = а 

Overlapping the plasmasphere is the ring current. 
Tt consists of a hot, thin plasma, with densities norms 5 
less than 18° m`? and temperatures greater than 19 keV. 
Again, явне is the dominant Species (935 of the er N 
density) while O* and He* are typically minor constituents: 


During quiet times the principal source or these и 


appears to be the solar wind. During active times, a large 
0” constituent indikates onespheriesssoureeesn {Rese 
6,221: | 

The plasmasheet is the region outside ene 


plasmapause, and like the ring current, consists ortam 1 
diffuse plasma. Densities оп the order o£ 14” = 13> m a K 
temperatures near l0 keV are common. At quiet times; 1 
composition 15 90% Н” and 18% Ot and Het. АР ає« үес 
large O* ratios (>50%) are observed. This is a mixing zone, 
containing particles of both terrestrial and solar “or pm 
Large negative potentials are found on spacecraft in the 


plasmasheet, particularly vhen the satellite is in eclipse. 


Cu 


Table 1 shows the, various ions in  theé* magnekosphere m 
thelr sources. {Rees= 6. лаве 
An interesting feature of these areas is tie 


Variability of the ion composition of the plasmas observed 


ви сте periods. In the ring current, Eor 
example, it has been observed that O* contributes a minimum 
of 23* of the storm time energy density, and may contribute 
much more. Also, the Не“ contribution increases, and 
adiabatic acceleration of these solar vind particles may 
leave them with energies up to 32 keV. Table 2 gives some 
КЕН] values £or storm time ring current Composition. 
ІКегс. 7,8] 


The plasmasheet region has a direct source of 07, іп 


ла aligned beams streaming from the ionosphere. The 
HOT ratio of the beams has been found to be 6/1  ІКеє. 
B. een e “has been observed that in the L = 5 - 


won e.g., the SCATHA orbit), during periods of high 
magnetic activity (K,*5), the O* density becomes comparable 
ww that of H*. (The McIlwain 'L' parameter is defined by 
R=L cos28, where R is the distance in earth radii, and 8 is 
ww magnetic а: де.) Ко. is ameasure of the general 
level of magnetic activity caused by the solar wind. The Ko 


И із quasi -lögarithmie, and runs from 9-9, with Ø the 


mew end and 9 the high. several occasions, the 07 
density has been found to reach 79%. ШаЕТоГ > а ови = оте 
EE > for the storm time magnetospheric bulk plasma 


ОШ Оз Стоп. (Refs. 6,7,8) 
3. Instruments 
The data for this experiment came primarily from the 


versity of California, San Priego auroral particles 


experiment,  SCOSE This experiment was Composed ot ap 


detectors}; two pairs of rotating ion and electron detece TEREE 


and one fixed ion detector. Figure 4 depicts the detector 
arrai ement: The rotating detectors 5саппед 11 оспе ЕН 
planes and vere designated the Hi set and Lo set. Тре Бе 


detector rotated from -20 deg to 200 deg, with 70 аед 2. 
the spacecraft spin axis. The Hi detector scanned the same 
range, but it vas more symmetrical,  vith its midpormn mM 
line parallel to бле р ийг S Thus, it could not depress 
more than 20 degrees belov the forvard end plane. The fixed 
detector looked radially away from the spacecraft,in the 
same plane as the Lo detector. The detectors had ап anc 8 
resolution Of S@ey 7 де The geometry is illustrate mM 
Figure 5. {Ref. 4} 

The Hi set covered an energy spectrum from M 
81 keV. The Lo rotating set, and the fixed detector, 
scanned from 1 eV Бо 2 ew. For all detectors, a compere 
energy scan required 16 seconds and vas covered in а еви 
of 64 logarithmic steps. The energy resolution at each step 
was approximately 205. Additionally, the detectors could be 
ordered to dwell at fixed positions and/or energies.{Ref.4} 

The detectors меге composed of three parts; an 
electrostatic analyzer, an electuostatic focussing. Lema NN 
a spiraltron particle sensor. The analyzer provided the 
energy differentiation for the system through curved “pi 


and an applied voltage. The lens focussed the particles 
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eues ensor center Dy means of two grids, one at ground апа 
encsotherzst the analyzer potential. The sensor was rated 
to 197 counts per second. {Ref. 4} 
4. Ion Gun 

An ion gun was installed on the SCATHA spacecraft to 
wiestigate the efficiency of an ion emission system in 
modifying satellite potentials. In particular, it could be 
used to develop negative voltages on the vehicle. The 
experiment utilized xenon gas ionized by cathode discharge 
and accelerated by either a l or 2 kV potential drop. The 
beam current could be varied incrementally from 90.3 mA to 
DU mA. The package also included an electron source that 
EEUd be configured to neutralize tne particle beam or to 


provide an electron beam. (Ref. 9] 


SCC OTHER SATELLITES 

Although SCATHA was the primary source of data for this 
thesis, two other satellites provided significant 
штоста јоп. These satellites, ATS-6 and ISEE-1, will be 
ес У а:5сц55ед 1п the following paragraphs. 

I ATS-6 

ATS-6 vas a large three-axis stabilized spacecraft 

WE in geosynchronous orbit in 1974. Its major features 
are shown in Figure 6 . The large antenna vas significant, 
since it fostered differential charging by shadoving other 
spacecraft surfaces. The satellite had an environmental 


measurements package above the spacecraft dish. unus 
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Figure 6. ATS-6 Satellite fRef. 12) 
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ПЕ ВЕЛЕНО еӨ particle detector sets similar to 
ешозе found on 5САТНА. (Ref.2? 


2. 15ЕЕ 





The ISBE-1 sstellitezwas Tauncehed in 1977 as part of 
8 xece spacecrirt mission to investigate the magnetospheric 
plasma. with this in mind, the satellite was designed to 
minimize absolute and differential charging by using 
Паола materials and iow impedance connections to 
ground. Regardless, the satellite was observed to charge on 
Several occasions. The spacecraft was equipped to measure 
electric fields, plasma waves, and plasma composition. The 
n importanc of these, for our purposes, is the plasma 
веро Сі оп experiment wissen <conrincd” an pon mass 


spectrometer. {Ref. 19} 


D. SPACECRAFT CHARGING 

The theory o£ spacecraft charging will now be reviewed, 
Ема: СгекепЕ1а1 charging will play a significant које іп 
Ши discussions. Іп general, a probe immersed in a plasma 
will develop some potential relative to that medium. = 
wwe пад" potential is determined by the balance of the 
mw enLts incident upon the probe. The primary currents in 
а па the potential are; the ambient plasma electron and 
шах, the photoelectron emission due to sunlight, and 
nr хез due to secondary and backscattered electrons. 
ИШ е currents are schematically illustrated in Figure 7. 


ШЕ eft half of the figure depicts typical eclipse 
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шоогоо (negace satellite), and the right half the 
Ls present in sunlight. {Refs. 11,12,13} 


тесте, епеееадоптлапс current 15 that of the 


mio oelectrons which, although material dependent, has 
ШОО Сат Saturation current densities in the range 5 - 50 
HA/m?. This current density overshadows that of the ambient 


plasma, which is generally below 1 BA/m? for electrons, and 
ВА? for ions. The value of the secondary electron 
emission current is highly dependent on both the material 
and the electron temperature. It is this material 
dependence which allows differential charging, even in 
eclipse. A normal secondary emission curve will have a 
yield greater than one in some middle energy range 
(typically 50-500 eV), and yields less than one at higher 
and lower energies. The point in the high energy range 
where the yield decreases to one is called the cross-over 
ewe and is critical in determining the effect of the 
secondary emission current. The energy at which this point 
decurs, and the peak yield, varies with target composition. 
Low energy ambient distributions (Ts5 kev} impacting the 
Wee craft Will create large secondary ‘fluxes, И Є з СС 
greater than the ambient electron flux, and the spacecraft 
will charge positively. rohr energy electron distributions 
(T25 keV) will produce lower secondary emission yields, and 
mm@em vehicle will gain negative charge. Thus, the potential 


m nc Satellite is determined by the high energy electron 


Т 


flux. The threshold electron temperature" tor which Chagas 
may occur is typically Іп Епе 5-19 Kev SENE There is an 
additional citeria, the critical energy, which also must be 
met. This is an upper energy bound on the distribution 
function resulting from magnetospheric convection processes. 
ТЕ is also termed the Alfven boundary, and must exceed 15-29 
keV for charging to occur. (Refs En ГІНЕ 

А simplistic equation to represent the situation may be 
се и о 
ШЕ, Тас “ олиш e İL 5 Е85 


TRUS, during daylight, when the photoelectron current 


dominates, the net current is positive (toward Еше 
spacecraft), and the spacecraft will start to charge 
positively. As this happens, the flow otf ions СО 01 


electrons from, the spacecraft will be inhibited, “angie 
value of Inet will decrease. When the value reaches zero, 
the currents are balanced and the spacecraft will remain at 
that potential until the balance is disturbed. For eclipse 
events, Јрносе із zero, and the balance 15 effecti TEn 
between the ambient and secondary electrons. For с 
energy electron distribution, the potential is restrictor Eea 
low values (+2 to +5 V), while for high energies ye. 
potential may reach - 1220 -2ЙоКУ. [Reier 27225 

In addition to absolute charging, a spacecraft ш s 
charge differentially. That is, separate insulated portions 


ОЁ the spacecraft surface may charge 120) aiffereme 


Potentials. ТЕ two insulators are exposed to the same 
environment, but have леве ве те зесопдаку emission 
characteristics, then their ability to shed the excess 
charge in the form of a current will differ, and a potential 
will develop between them. шїї ээг ır Бо“ portions аге 
exposed to differing environments, they will each charge to 
ши cerent potentials. The environment may differ over the 
range of a spacecraft in a number Of ways. Frequently, 
pitch angle anisotropies in the ambient plasma may cause 
different fluxes to impact separate spacecraft regions. 
Additionally, if опе portion of the spacecraft shadows 
another continually (see ATS-6), the shadowed area will not 
emit photoelectrons and will charge negatively compared to 
the rest of the satellite. Also, (НЕ ончет “Orbits, ап 
electric field induced in the satellite frame of reference 
(E=vxB) may cause variations in the ambient particle 
пи ii butions. {Refs. 2,12,13} 

Differential charging is believed to be the cause of 
mumerous control logic upsets on various spacecraft. As the 
differential potential increases, there will be a tendency 
шээс to discharge іп the form of an arc. This arc creates 
an electromagnetic pulse which is read by the control 
circuits as a command signal. Borslargergdiscehargesz there 
may even be physical damage or destruction to the arcing 


areas. {Refs. 17,18} 
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Spacecraft charging may be controlled by active and 
passive measures. Active trials using ion and electron guns 
were carried out on ATS-6 and SCATHA with some success. 
Passive attempts, such as those on the ISEE satellites, 
using low impedance connections and conducting surfaces were 


also moderately successful. {Refs. 19,19,29} 


E. SPUTTERING 

The general theory of sputtering is very significant to 
this thesis, as it is postulated to Бе the source of the 
observed spacecraft generated ions. Sputtering is defined 
as the ejection of material from solid surfaces 0475 
bombardment. The observed features of the sputtering yield 
(number of target atoms sputtered per incident ion) are; а) 
its dependence on atomic Pr pem b) its correlation туш ћи = 
periodic table, с) its correlation with the <a 
sublimation of solids, dd) its relation to the еї ТЄТЄ M 
momentum transfer, e) and its dependence on incident ion 
energy, such that it increased to a maximum and then 
decreased. Figure 8 isa plot of the sputtering yem 


aluminum as a function of incident ion energy and mass (hee. 


21]. We see that the yield maximum increases with ion mass, 
and shifts to higher energies. These observations led to 
theories based on atomic collision cascades. Such cascumes 


are created by one incident ion impinging on the lattice jana 
transferring energy со other 3505” which in turn 


рагсїстра ё” сое. These theories are based on the 
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Figure 8. Sputtering Yield for Aluminum {Ref. 21} 
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simple elastic collisions of clai cal mechanics and 
therefore satisfy the observed dependence on atomic number 
and momentum transfer. The correlation to the heat of 
sublimation is explained by assuming the effect o 210 
surface binding energy. The maximum in the energi 
dependence occurs by postulatini some Tape men energy 
deposition depth, beyond which the ability of the cascade 
particles to reach the surface decreases. {Refs. 21,22,23} 
The most comprehensive, detailed, and successful of the 
collision cascade theories is the Sigmund - Thompson 20100 4 
This theory is actually two theories developed separately, 
but they complement one another by explaining а 2 2 
facets of sputtering. The portion developed Бу 5190 
treats the sputtering yield of ал amorphous targe n 
Thompson's theory explains the energy distribution e W 
sputtered atoms. The development of both these theories 
will be briefly discussed, followed Бу а comparison with 
experimental results. {Refs. 22,23,24,25} 
1. Тре 5тлатила Theory 
The Sigmund theory focusses on the amount of ey 0 
deposited in the surface layers of the solid as the driving 
factor Боп eme e og The basic expression ff co НЫ 
sputtering yield із then S=BF, where» г is a Соён Б 
depending only on the properties and state of the target, 
and F is a function describing the interaction of (Ле 


particles. In order to evaluate this expression, Sigmund 
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СӘС as 1:15 Starting Point the sputtering of an amorphous 
Eust) cand used transport theory to describe the collision 
cascade in random media. (Refs. 23,24} 

The theory assumes a planar surface and an atom 
starting its motion at time t = 0 and position x = Ø. The 
ИЗТ quantity of interest is G(x,%,V,t)d’vo dx. Bons 
is the average number of atoms moving at time t іп а layer 
EN with avelocity (Vo,d’veo) due to an impact by a 
Ele atom at velocity У. This guantity is then used ina 
Boltzmann transport eguation that eguates the initial number 
G with a later G which is dependent on the probability of 
1 1510п. The equation is then integrated to eliminate 
some of the variables. Two functions are defined to 


= 1 1 Су the integration. 


>) F(x,Vo,V) - ca dd 
(1-3) Mv = | My аза эа 


The function F is the total number of atoms which penetrate 
EN але х місп а velocity (Vs,@2ve) during the collision 
cascade development. H(x,V) is the sputtering yield at the 
plane x for a source of sputtered particles at x - 8. Mae 
function may be defined for either forward or backward 
epn tering. For our purposes, we will choose backward 
EEECLtering. Бае су сіні весла 15 emission directed 


eepos ite to the initial direction of incidence.íRefs. 23,24) 
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The remaining expression in H is expanded in terms 
of Legendre polynomials and then transformed into a moment 
equation by multiplying each term by x°” and integrating over 
ке With the assumption that the electronic losses to 
electrons are small and isolated Екот nuclear collisions, 
the electronic and collisional cross sections may be 
separated {Ref.26}. This yields a final expression which 


includes the electronic stopping cross section, and the 


differential cross section for elastic scattering. xu 
22,24) 

The next step “іс to determine the function 
describing the electronic sand elasScici oss ccc The 
electronic cross section used is Se = K БЕ“, where = 
constant K is determined using Thomas-Fermi arguments 
(нео ти The elastic cross section is derived from 


classical mechanics, and a power approximation m 


Thomas-Fermi  interatomic potential. The equation is then 
given Dy: 
22-12 ОСТ) - СӘН == от 


where C is a constant, T is the transferred energy, and E 
the initial energy. The variable m comes from the power 
approximation. m-1 corresponds to Rutherford! scattering mmm 
m=@ approximates scattering from a Born-Mayer potential. ТЕ 
these values are substituted in the expression from the 


previous paragraph, the following is achieved; 


a) ЕСХ,ШН Р) - 3 ЕХУЕВ)- 

412 М Со Оо 
where ß is the direction cosine of the ejection vector, and 
ELE B) is a function describing the distribution of the 
energy deposited in the solid. Co is a constant, and Ue is 
the binding energy. bem erer ye che tunctlon F(x,E,p) 15 


found to obey the following relation; 
1-2 Е(х,Е,8) = «a N Sn(E) 


where « is a dimensionless constant dependent on the ratio 
М-/М:, N is the atom density, and S„(E) is the nuclear 
stopping power, given by Sa(E) = | PEJO. Additionally, we 
шэг restrict ourselves to sputtering in the unit normal 
Шігесбіоп. The function. «= is illustrated “in Figure 
ew Refs.23,24) 

Using all of the previous relations, we can then 
define some general expressions for the sputtering yield. 
Sere low energies (EslkeV), we have 
(1-7) Б(Е,В) = 3 « 34MiMz E_ 

Ап 2 (MitMa)? Ue 
Sa(E) in this equation was evaluated by integrating do with 
a value of m=@. The Born-Mayer potential is more accurate 
at lower energies than the Thomas-Fermi approximation. The 
гс топ tor higher energies is; 


ı ) SCE, B) = 3.56 Q Z М, SE | 
er) (Мі Мо) Ја 


279 





Figure 9. «(Mi/Mi) уз. М,/М, (Ref. 23] 
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m wtunctlon 5. (Е) го спе reduced stopping power. This isa 
universal function for  Thomas-Fermi interactions (Refs. 
24,27). The variable e is the reduced energy, defined by 
€ = aM;E , where 7; апа Zz are the atomic 


(ZiZ2e2 (MitMe2) ) 
number, and a is the Lindhard screening radius. Tabulated 





values for s„(E) are shown in Table 4. The sputtering yield 
dependence on В vas roughly determined to be  (cos(8))-*. 
This dependence was derived  empirically by Бата. 
Moers .22,23,24,27} 
2. The Thompson Theory 

The Thompson Theory was a complete theory of 
sputtering published in 1968. The portion involving the 
determination of the sputtering yield was overshadowed by 
the Sigmund Theory, but the energy spectrum determination 
Mas endured. Thompson started his derivation by assuming an 
Infinite solid with recoils at energy Ez from some primary 
event. Then, with a density of recoil atoms q(Ez)dEz per 
unit time, each generating their own collision cascade, the 


total number of atoms slowing down through some energy range 


ЕМ 15 
2-9) n(E') = Г cq UEM 2 ава 

E 
The mean rate of energy loss is given by VdE'/dx, where V i 
Ww <evelocity. Assuming an isotropic distribution, опе may 
define 


E I0) n(E',rZ)dE'd&'-f* q(Ez)Uu(Ez,E') dEs dE' dga' 
= У dE'/dx 4 т 
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From this, the flux may be determined, 
RUNE Ф(Е",Г)4Е"40"- 9 п(Е!,Ү)созӨ4Е"4аа! 


The theory now assumes that the infinite solidas 
cut in half and the £flüux G is observed Also, some 
expressions are needed for the functions above. U(E2,E') is 
the number of displaced atoms and may be approximated by 
BE2/E', where ñ is a constant of order unity. Similarly, 


dE'/dx may be replaced with E'/D, where D is the interatomic 


Spacing. Using these relations we obtain 
(1-12) Ф(Е",С)4Е"48"-3 Ор Гэа(Е-)Е-аЕ» созе да" дЕ! 
E'? ы! 4 к 


This function éó is the flux inside the surface, and must be 
transformed to that observed outside. This is accomplished 
by assuming a binding force normal to the surface which will 
decrease the normal kinetic energy, but leave the parallel 
component unaffected. Thus, the effect will be a bending of 
спе particles trajecrer With the energy equations, and 
the relationship between the inner and outer angles in hand, 
we come to the final expression 

ÇE) Ф(Е,Ф)дЕдо= В Р cos |” а(Ег)ЈЕгдЕљ адЕда 

Ат (1+Ј;/Е)ЗЕ2Ј 

In this eguation, Us is the binding energy, E is the emitted 
particle energy, and Æ is the angle of the trajectory from 
the surface normal. The chief features are the 


proportionality Com ЕЛ ОВ = дерепдепсе. 


2:8 


Целе 10 15 з рјоЕ о this function versus energy. For 
high energies, the spectrum will fall off as 1/E2, and for 
low energies, the binding energy will have an impact. The 
maximum of this curve occurs at Ue/2. Figure ll isa plot 
Or the measured energy distribution,with the model overlaid. 
It can be seen that up to approximately 1000 eV, the model 
matches the observed spectrum quite closely. Since the 
majority of the sputtered particles are in the energy region 
that follows the l1/E? curve, this theory is quite effective 
in describing the data. {Refs. 22,24} 
3. Comparison of Theory with Experiment 

To evaluate the applicability of these theories, we 
must determine how closely they recreate experimental 
results. Therefore, we will look at the sputtering yield as 
ENEUnction of ion energy, angle of incidence, incident and 
target masses, and surface binding energy. The energy 
dependence of the sputtering yield is contained in the 
energy dependence of the stopping power. Calculated yields 
using equations 1-7 and 1-8 are compared with experimental 
os in Figure 12 for copper (Ref. 23}. The solid line 
is the high energy result, and the dotted line is the plot 
for the low energy equation. The experimental results are 
marked Бу the various geometric shapes. п addition 
calculated and measured yields for various materials are 
ШО ра іп Table 5 {Ref. 23}. It can be seen that the 


correlation is reasonable, normally within a factor of two. 
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Some of the experimental yields are questionable, as the 
sputtering yield is known to have a dose dependence. The 
Sigmund Theory also tends to overestimate sputtering at very 
low energies («1200 eV). In most cases, experimental data 
indicates that the yield goes to zero in the 54-100 eV 
region.{Refs. 22,23,24} 

The sputtering yield is a function of Mi aik CE 
and in the factor a(M=/Mı). Due to the above mentioned dose 
effects, yield data for masses are generally normalized to 
the self sputtering yield or the argon sputtering yield. 
Data for various targets are presented in Figure 13 {Ref. 
2318 The top plot displays results for impacts on sim 
The mass of the sputtering ion in AMU is along the 
horizontal axis and the yield of the ions normalized to the 
argon sputtering yield on the vertical axis. The other 
plots are similar for copper and gold. The agreement with 
theory for silicon is quite good, but decreases "POEM 
heavier targets. Sputtering yields Ёок Н" апа4 о“ оп ске 
are illustrated in Figure 14. Since in general we will be 
discussing silicon or silicon dioxide as a target, wen 
disregard the lack of agreement for heavier targets in the 
model which we develop. {Refs. 22,23,24} 

The angular dependence of the sputtering yield is 
Considered to ye propor: oni Seo ees on ire The exponent, 
n, is a variable quantity ranging from a little less than 1 


to 1.7. For the situations we will be discussing 580 
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Figure 13. Mass Dependence of the Sputtering Yield 
for Various Atoms {Ref. 23} 
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M=/Mı > 1, the value of n is generally taken to Бе ШАН 


Figure 15 15 a plot of the dependence on angular incidence 


from experimental data (Ref. 23}. The dependence ranges up 
to an incident angle of 70 degrees. Beyond this poin r p s 
sputtering yield decreases rapidly to #0. If we attempt to 


model using only the 0-70 degree range, we will probably be 
underestimating the yield by approximately 20%. {Refs. 
22 99 24] 

The yield depends directly on the binding energy, 
thus the determination of this energy is quite important. 
Sigmund suggested that the appropriate value for the binding 
energy would be the sublimation energy of the material. 
This generally leads to sputtering yield values that are 
slightly аа. Recent work has given theoretical backing 
for a higher value, but no definite numerical quantity was 
put forward {Ref. 28}. It has been suggested that the value 
providing the best results is 17/10 of the sublimation 
energy íRef. 29}. If this value is used to recalculate the 
numbers in Table 5, we qet the results in Table 6. Ї гэ 
obvious that in general, the calculated yields in Table 6 


are much closer to the experimental yields, than those 


calculated for Table 5. For calculations in this thesm mm = 
will use the larger value of the binding energy. ЕЙ 
О! 


=> 5(8) / 500) 


Е = 1.05 кем 





Angle of Incidence 


Figure 15. 
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The Dependence of the Sputtering Yield 
on Angle of Incidence {Ref.23} 
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4. Limitations oM the Sigmund есу 


The theory was devised for amorphous monatomic 
targets  bombarded by atomic ions. It is not so reus 
applicable to the bombardment of molecular solids, as will 
be necessary in later calculations. Since several of the 
reference data points are concerned with bombardment with 
molecular oxygen, it will be relevant to discuss this type 
of impact. The case of molecular ions on molecular targets 
may be treated through computer simulation, but it is an 
extensive, complicated project, as the probabilities for 
each collision combination must be calculated. In most 
cases, it is more appropriate to utilize experimental data 
if it is available. Such data is often difficulter 
obtain. Only recently has an effort begun to systematically 
tabulate the known experimental results and compare them 
with the theory. The presently tabulated data concerns 
atomic ions on monatomic targets {Ref. 38}. A great deal of 
data is also distributed randomly through the literature on 
other ion-target combinations. 

An additional bit of information is available for 
the yields of monatomic molecules. It has been founda 
several atoms that the sputtering yield per atom ЕСЕ 
molecular ions is slightly higher than that for atomic 0 
(Rer 311. The results are listed in Table 7. The yim umi: 
are greater by approximately a factor of two, (wi e ааа 


silicon they range fron по И, There is a 
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о кепсе. but Yt is пос drastic. While it is difficult to 
accurately model the more complicated situations, it should 
be possible to obtain some approximate results for monatomic 
molecules that are within a factor of two of the actual 


results. 


F. NASA CHARGING ANALYZER PROGRAM 

The NASA Charging Analyzer Program (NASCAP) was 
developed to accurately model the dynamics of spacecraft 
response to realistic plasma environments. NASCAP can 
effectively simulate the charging, in both laboratory and 
magnetosphere environments of objects ENaC are 
Geometrically, materially, and electrically complex. The 
program utilizes a timestep procedure, calculating quasi- 
static steps. The dynamics are driven by charge 
accumulation from external sources, charge depletion, and 
Semauction in dielectrics. Each timestep includes a full 
three dimensional electrostatic potential calculation. The 
calculation of Poisson's equation over a grid is time- 
staggered with a procedure in which incident charged 
particle fluxes, leakage currents, emission currents, and 
induced spacecharge effects are computed based on the 
Pure Quasistatic conditions. NASCAP is limited to 
Situations where the Debye length is long compared to the 
dimensions of the examined object. This COnd1elon is met 


moms the SCATHA spacecraft. NASCAP can output a time history 


of spacecraft charging, potential contours, charge come | 


current contours, and particle trajectories. 
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II. OBSERVATIONS 


А. PREVIOUS OBSERVATIONS 

шапа taken from three satellites,  АТ5-5, ATS-6, апа 
ПЕ - 1, mili be discussed, rior to presenting Р78-2 
(SCATHA) results. Observations of the three primary 
phenomena (triangle peak, shadow peak, diffuse background) 
on satellites prior to SCATHA will be illustrated. 

ШОК АТ5-5 

The presence of ions below the charging peak was 

ww reported for ATS-5 {Ref.l}. АТ5-5 data taken in 
eclipse in 1969 and 1979 showed ion fluxes ina broad range 
of energies below the charging peak. A typical example of 
such data is shown in Figure 16. The data is presented in 
grey scale spectrogram format, a common method of displaying 
the data concisely. The figure is vertically separated into 
three components. The top primarily displays magnetic field 
information, unnecessary for our purposes. The remaining 
two sections indicate the count rates of electrons and ions 
in the detector versus time and energy. The horizontal axis 


defines increasing time, and in this case is labelled with 


ШЕН nours of day 274 of 1970. The verticai axis 15 
associated with the energy of the particles. The minimum 
(50 eV) of the energy axis is located between the two 


sections, and energy increases upward for electrons, and 
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downward for ions to a maximum of 50 kevV. The« magnitude of 
the detected countrates is indicated by the shading of the 
Ш сс from black to white, with black corresponding to zero 
eounts, and white with maximum. l pene ~adetector” is 
 Ппғагей at а certain time in an energy channel, che 
spectrogram will contain a black spot at that time and 
energy. A feature such as this can usually be distinguished 
from a zero count rate because saturation points are 
normally surrounded by intense count regions, revealed as 
bright white areas. 

The data displayed in Figure 16 are from day 274 of 
Шола. The satellite enters a region of hot plasma at 9450, 
as indicated Бу the increased intensity of high energy 
(24500 eV) electrons. The satellite charging peak may be 
seen as the intense white region in the proton spectrum 
between 0620 and 0728. From the energy of this region, we 
can infer a satellite potential of approximately .75-2 kV. 
The black spots on the peaks are an example of how detector 
saturation is displayed. Шеноа Бес киш oF” ous is 
evident in the energy bands below the charging peak as the 
grey area that clearly extends from 1 kV to about 399 V, and 
may even continue to lesser energies. Шесе particles dic 
not appear during every charging event, so the source was 
difficult to determine and no complete explanation was put 


mor Ward. 


2. ATS=5 
ATS-6 also showed evidence of non-ambient source 
contamination of the ion data. Secondary or shadow peaks 


were observed at energies belov the inferred satellite 


potential. An example of such data is illustrated in Figure 
ee This figure is a grey scale spectrogram similar to the 
one shown im Pigure 16. The data presented are from the 


North-South detectorzon АТ -5 сп зу 2458 сс ован 

This plot displays two phenomena. The first is the 
ion peaks, which are clearly presented elsewhere {Ref. 2}. 
These ions were denoted "spacecraft generated ions". The 
second is an apparently analagous phenomenon in the electron 
spectrum. These observations are the so-called "Minnesota 
spors”. In addition to the electron and ion spectra, this 
figure displays the detector pitch angle in the top section, 
which is related to the rotation angle. 

The spots are seen in the electron spectrum between 
1040 and 3080 eV as intense white points that appear to form 
triangle shapes. If one compares the periodicity oe 
electron triangles with the detector look angle periodicity, 
a direct correspondence is evident. As tne top of the 
spacecraft, where the detectors are located, was in shadow 
at this time, the source electrons were determined to һе 
secondaries from the surface of athe vehicles These 
electrons had been emitted from a surface differentially 


Charged with respect to the detectors, and thus the 
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Figure 17. ATS-6, Spectrogram, Day 248 of 1974 {Ref.2} 
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charged with respect to the deteetors, and thus the 
distribution had been accelerated by a potential асори 
The source of the flux was eventually determined to be 
secondary emissions from the University zei Minnesota 
detectors also located on the satellite measurements box. 
IRef. 2) 


б) SR 





Negative charging events on the ТЗЕЕ за 21. 1 сав те о 
relatively rare, as the satellites were specifically 
designed to avoid the problem of spacecraft charging zı 
1 cid, however, charge to significant negative levels оп а 
few occasions. During one such event, analysis of the ion 
mass composition revealed some interesting features. Figure 
18 is a presentation of the ISEE-1 mass spectrometer reading 
from March 17, 1978. The vertical axis is a count za: ше 
the horizontal axis is the applicable mass channel. The 
mass channel is related to the atomic mass by, 


(251) AMU = 2010 1 
MasseChanmne lw. — 


Thus, mass channel 54 corresponds to H*, 15 to N*, amda M 
Oetker.) 
АС 0633 UT (Universa? Time), the instrument 


responded with a broad peak in the high Atomic Mass Unie 


(AMU-30-100) channels. This peak had a maximum in channel 
L0 Corresponding toe > А: There was also a peak at 
channel mi, п 56 The normal value of the ambient atomic 
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Pigure 18. 


oxygen is indicated by the solid black= line under Chem ни 
The data from 9721 UT are quite similar, except that saw 
point the maximum has shifted to channel 11, 28 AMU. Again, 
this is larger than the indicated ambient oxygen. At the 
time these observations were first made it was concluded 
that these molecules were probably hydrocarbons and other 
QUEGaSS па зеоан не 5. In light of spacecraft surface 
sputtering, however, опе may postulate a different flux 
source. That is, these particles are ions emitted from the 
spacecraft surface. Since the surface is primarily silicon 
dioxide (S103), the ejected particles would largely Dom iE 
O, and Oez. The masses of these are 28, 16, апа 32 7740 
respectively, numbers which would fit nicely in = 


observed data. {Ref.19} 


ВТМ 13-22 SCATHA) 

Observations of spacecraft generated ions are presented 
for five days. The first two are for eclipse charging 
events in 1979, the second pair are for ion gun 1 08 ёё 
charging, апа the last is for eclipse charging in и 


Jo Day on Э 


This first eclipse event illustrates the triangle 


ка а 
{N 


peak and its associated angular dependence. Figure =i 


2 


the grey scale spectrogram for a one hour period of G m 
OE T979: This is a composite of the electrons from eşe ni 
Getector, and the ions from all three detectors. bach zo: 


the ion displays has zero energy at the top and increases 


energy downward. The electron presentation is reversed, 
with low energies at the bottom, increasing upward. Between 
17749 and 1820, the triangle structures can be seen in the Lo 
detector data. From the Hi detector data, it can be seen 
that the spacecraft charges negatively to about 244 у at 
1750, decreases to 100 V at 1815, and increases rapidly to 
EJ V at 1825. Multiple shadow peaks are visible during 
this period in both the Hi and Fix detector data. 

Ere ole: or the “log. ton distribution 
Ennction (phase “space density) versus energy for the Hi 
detector. Dnemsoseecram- is charged to -—259 57, in eclipse, 
as indicated by the peak in phase space density (f) at 258 
ev. The shadow peak is the secondary peak at 78 eV. At 
enis time, the detector is parked at 272 degrees, 
approximately parallel to the spacecraft spin axis. 

Figure 21 is a diagram of the spacecraft potential, 
triangle peak energy, and Lo detector look angle versus 
time. The energy of the peak varies directly with the look 
angle, with the energy minimum occurring at the angle 
maximum, that is, looking down and away from the spacecraft 
(where it approaches closest to the spacecraft surface). 
The shadow peak energy maintains a relatively constant ratio 
with satellite potential. In this Case, Ep/Et is between 6 
апа 7, where Ер 15 the satellite potential, and Et is the 


triangle peak minimum energy. Gn Sns significant, аз ме 
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Peak Energy, 


= 2 


would expect a differentially charged surface to maintain a 
constant ratio over time with the satellite mainframe. 
2. Bay 96 ОБ 19173 

The second eclipse charging event, day 86, гэ 
ма Јах исо that of day 83. Figure 22 is a grey scale 
spectrogram of the ion data for a one hour period on day 36 
EE 1979. The triangle structures are again visible from 
1645 to 1749 in the Lo detector display. From the Hi data 
we see that charging begins at approximately 16540, апа the 
potential varies widely throughout the ое гос. А 
fluctuation in the satellite potential due to a change in 
the electron Alfven boundary is seen at 1705. Tha 
Spacecraft achieves a maximum potential of -600 Vat 1725. 
During the period, a shadow peak is clearly visible at 
energies lower than the charging peak in both the Hi and Fix 
detector data. 


К лс з К lots the ion dis соц Erom the Hi 


detector. The satellite potential is -34@ V as indicated by 
the charging peak. There is a well defined shadow peak at 
1@@ eV. In addition, there are various other small peaks 


oelow the charging peak that may indicate other shadow 


peaks. Frequently, other shadow peaks appear in tne 
EpGctrograms, maintaining a constant separation с 
between each other, and varying as the main peak varies. 


The detector is again parked parallel to the spacecraft spin 


exis. 
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The strong dependence between the triangle pem 
minimum energy and spacecraft potential is 111и5бгтаСеа 0520 
Figure 24 for this more variable day. It 1s clear that the 


minimum energy of the triangle peak increases with the 


increasingly negative vehicle potential. The ratio, крш ве 
is maintained between 6 and 6.5. Again, the energy of the 
Shadow peak varies directly with the detectcr angle, 


achieving minimum energy at the detectors maximum angular 
extent. The shadow peaks usually lose their definition 
below 112 degrees. 


These two examples illustrate typical обзекуа и 


dur ing eclipse passages for which negative charg ine 
occurred: ву comparison? Wich 4725-65 we would expect 
similar observations during daylight negative charging 
events. ATS-6, however, was stabilized and had its large 


solar array to provide a shadowed area which would simulate 
an eclipse. SCATHA had no such array, апа in addicromm = 
spinning so that no area was shadowed long enough to charge 
to a high level. Observations of daylight charging | 
that the spacecraft potential was highly spin тоа ни 
making it difficult to use these data for this We a 
However, induced charging events with the ion gun prov 1 
additional тападағог даауітанһі азса” Two examples  ocr БЕ 


experiments are presented. 
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3, Day ЭА оинаш 


The first example of ion gun induced charging ч: 2 
94 ge 121727 The ion gun was in operation at I KV ame ГУ. 
Figure 25 is the grey scale spectrogram for a 1 hour р o s 
from снас “ға The data presented are from two ion 
detectors, with the Lo detector in the upper diagram, апа 
the Fix in the lower. The ion gun is turned gene s 
approximately 1405 and the spacecraft charges promptly to a 
final value of -60 V. Shadow peaks occur at ІІ and 20 eV in 
the Lo detector between 1420 and 1430, and between 1440 and 
1455. 

A curious feature is the appearance of twin shadow 
peaks in the Lo data at 1444, with bands between in which no 


ions were detected. When this phenomena began, the detector 


ba 


had just been parked at an angle of 28 degrees to th 


iD 


Spacecraft spin axis, away from the satellite body. on 
detector had previously been scanning. When the data 
printouts are examined, three bands clearly emerge; one at 
11 eV, another at 22 eV, and the last at 38 eV. We may 
infer that certain energies, or trajectories, do net reach 
the detector at that angle. These events indicate something 
about the location of the emitting surfaces. From the 
striations in the Lo data, we would expect non-emitting or 
non-differentially charged surfaces to be sandwiched between 
emitting regions, thus giving the dark regions where по m 


are detected. Importantly, the Hi detector is viewing in a 
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Figure 25. SCATHA, Spectrogram, Day 94 of 1979 
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similar direction е and gives a similar 
response. The Fix detector, viewing 70 degrees from the Lo 
and Hi, reports a completely different environment. ІВ 
does, however, show a similarity in that a large energy gap 
exists. This indicates mixed differentially and non- 
differentially charged surfaces as above. 
4. Day 28098 of 1979 

Our fourth example is during an ion gun operaci TE 
1 КУ, 1 тА. Figure 26 is a grey scale spectrogram for a two 
hour timeframe. The data are taken from the Hi detector. 
In this figure, the spacecraft is charged to approximately - 
700 V, except during the two five minute periods when the 
ion gun was in "trickle" mode. Trickle mode is defined as 
the gun in operation with the accelerating voltage off, but 
currents 2 арр11е4 to the diodef Thus, the generated 
ions leak out at low energies. Triangle peaks are clearly 
visiblemrmeries'on display. This day is one of the few in 
which these peaks become evident in the Hi detector. The 
ratio, Ер/Ес, for these figures is in the 9-10 “аа 
indicating that the Hi detector is looking at an emission 
point different from that seen by the Lo detector on days 83 
and 86 (or different processes are at work). 

Figure 27 shows the distribution functions when the 
Lo and Fixed detectors are looking approximately parallel, 
radially away from the spacecraft. The Hi detector is near 


its maximum rotational angle. The satellite potential is 
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Figure 27. Day 200 of 1979, Ion Distribution Functions 
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Log f (S97 km? ) 


approximately  - 708 У. The Hi and Fix detectors exhibit a 
broad array of low energy ions, but such a spectrum is not 
evident in the Lo detector. The Lo detector does show a 
sharply defined shadow peak. The data have a similar form 
wreugnout this period, with subtle variations in the 
amplitude and energy of the diffuse spectrum, particularly 
with angle. 
о Бам 92 o£ 1981 

Figure 28 is the grey scale spectrogram from a 1981 
eclipse charging event. The figure displays 4 hours of data 
meome the Hi detector for both ions and electrons. The 
eclipse charging event is visible in the illustration 
between 0810 апа 8945 UT, with potentials varying from -1900 
Со -400 У. During this period, a definite shadow peak is 
also evident, in the -68 to -158 V region. 

Figure 29 is a Por e LOnmecn ser woution 
we tLons for day 92. me ова ве ве 5 :оп function has a 
charging peak at 0820 ас 90 еу with a value of 6.97. TRIS 
value is comparable to those on days 83 and 86 of 1979. 
Thus the detected countrates are as high as those 
encountered early in the mission. This indicates that the 
source of the detected particles is not dependent on 
satellite lifetime, and hence that outgassing is not a major 
eomer lbutor to the observed satellite generated ion 


phenomena. 
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Lit. DISCUSSION 


A. DPOSSUBEEOPONESSURGCES 

The ion source must be able to explain three different 
phenomena, ће triangle peaks, the shadow or mirror peak, 
and the diffuse spectrum. The prime feature that the 
observations have in common is the appearance of the subject 
ions at energies below that provided by the spacecraft 
potential. Thus, they must be generated in regions near the 
spacecraft where they do not experience the full  plasma- 
satellite potential. There are two prime candidates for the 
source, outgassing and sputtering. Outgassing is the 
emission from the satellite of contaminant molecules and 
atoms with energies near the temperature of the satellite. 
These are primarily trapped hydrocarbons and other 
atmospheric molecules. The outgassing rate may be increased 
by the inclusion of efflux from the attitude controma ex 
and the ion beam system. Sputtering is the emission of 
charged and neutral atoms from the satellite surface due to 
ion bombardment, and its theory has been described 
previously. A qualitative discussion of both possibilities 
follows, followed by a numerical calculation of the number 
of ions that would reach the detectors due to Sigmund- 


Thompson sputtering from differentially charged surfaces. 
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1. Qutgassin 
As the contaminant particles leave the spacecraft 
Pew are subgect to ionization from photons, and collisions 
with ambient ions and electrons. As the molecules. are 
ionized at various distances from the spacecraft they will 
return ina broad range of energies. This is а possible 
explanation of the observed diffuse spectrum. The question 
is then, how many particles outgas, and what percentage of 
these are ionized? The problem of outgassing and the 
outgassing yield has been studied previously fRefs. 32,33}. 


The dominant mechanism for the ionization of the outgassing 


molecules has been determined іп sunlight to be 
peotoionization, and in eclipse to be collision with 
electrons. The electron ionization rate is approximately 


two orders of magnitude less than the photoionization rate. 
In eclipse, however, the photolonization rate goes to zero 
and the electron collisions become dominant. In all 
probability, the electron ionization rate will increase 
slightly in eclipse, as the electron density will increase 
due to space charge effects. The ionization rate, however, 
will not exceed that of the photoionization. (Refs. 2,32) 
The distance within which the ions may be returned 
to the spacecraft will be the width of the surrounding 
plasma sheath. The sheath distance at geosynchronous orbit 
is approximately 100 m. With this information and the 


ionization rates above, it may be shown that, at a maximum, 
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approximately .02*$ of outgassing silicon atoms VAA 
ionized. In a central force field, with сопвегсуаста 55 
angular momentum and energy, calculations reveal that all 
the particles ionized in the sheath vill be returned 60 the 
spacecraft. The final result then, is that .92% of the 
outgassing molecules will be reattracted to the satellite. 
[Ref. 192 

When the spacecraft is first launched, the  ionized 
outgassing flux may be sufficient to provide a Sigma s m 
number of particles at the detectors. Outgassing, however, 
is an extremely time dependent phenomenon. Once the 
particles leave there is no way to replenish them, so the 
outgassing rate must decrease with time. Various 
experimental data have placed the outgassing as proportional 
to t^^, where п varies from 0.8 to 1.6, or as prope mi 
to exp(-t/41.6), with t іп days {Refs. 32,33}. Те ши 
e-folding time then would be 41.6 days. Even this would 
lead to a decrease in the outgassing flux by three orders of 
magnitude over the period of a year. As noted in day 92 of 
1981 above, however, the shadow peaks have not diminished 
over time in any significant way. This leads to the 
conclusion that the outgassing flux may be a contributor, 
but is not the prime source of the ions. 

What effect do the thrusters and ion gun have on 
this conclusion? An experiment was performed aboard the 


SCATHA spacecraft to measure the contamination (impurities 
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mixed with the hydrazine fuel) on the satellite due to the 
thruster operation. Two sensors were placed on the 
satellite, one in the bellyband, and one on the forward end. 
The sensors were capable of detecting 5 ng/cm? particle 
deposition and a current as low as 197+? A. The thrusters 
were located on the aft end of the vehicle, and were fired 
once per week during the bulk of the flight time to correct 
for precession. The results of the experiment indicated 
that no measurable flux of thruster contaminants was 
returned to the spacecraft. {Ref.34} 

The ion" gun is also ап ми Грке by source Or 
contaminants to the detectors. As noted above, it is 
located on the aft end of the spacecraft and directed away 
om it. During most normal modes of operation, the energy 
and angular momentum of the exiting particles would prohibit 
them from reaching the detectors. There is evidence that 
some ions are returned to the detectors, but they are 
observed at energies greater than the satellite potential 
(Ref. 35}. There is a possibility that ions emitted during 
trickle mode, at low energies, may be reattracted, and 
Wwe iribute Со the flux observed at the detectors. This is 
pusbably only a small intermittent component, as the 
operation of the ion gun in the above mode was not 
continuous or frequent. Thus, cae soUcgassing ~cElux is 


primarily composed of the emitted contaminants. 
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With just the emissions of the contaminant tla rt 
is still numerically conceivable that the outgassing flux 
could produce sufficient particles to account Wr UE 
observations. However, it is certain that the outgassing 
flux does decrease over time, and this is not observed in 
the data. Additionally, it is hard to conceive that the 
Satellite potenials are such that the randomly ionized flux 
is focussed into a beam one energy channel wide at the 
detector. Moreover, it is unreasonable that the focussing 
should be such that a small change in the detector look 
angle yields such a definite change in the observed energy, 
as in the phenomena of the triangle peaks. In sum then, it 
does not appear likely that ionization of outgas products is 
the source of the observed flux. 

2. Sputtering 

When this work began, it was initially assumed that 
the ion source was outgassing. As the work progressed, it 
became apparent that there were serious drawbacks in this 
assumption, as noted above. while searching for 
alternatives, the advantages of sputtering as a source 
became clear. First, besides a small dose dependence, there 
is no reason for the sputtering yield to decay over time. 


This attribute was very important in light ОЁ the 


significant shadow peaks observed in 1981. Second cme 
emission of sputtered particles from a differentially 
charged surface would lead to observed beams at the 
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detectors. As the particles left the surface with their 
various initial angles and energies, they would be 
accelerated by the existing electric field into an almost 
parabolic orbit. Thus, for each energy, only particles 
following specific trajectories would reach the detector, 
limiting observation angles. Since the detector energy 
channels have a resolution of 20 5, the bandwidth of the 
channels increases with energy. Therefore, a peak which 
would appear broad when distributed over the lower energy 
channels, can be compressed into a narrow peak in the high 
channels. The low and high energy channels of the SCATHA 
detectors are listed in Table 8. The similarity between the 
observed data and the data from the Minnesota spots also 
became evident. It seemed reasonable that if such a process 
ша оссиџе for secondary electrons, it could occur for what 
were essentially "secondary" ions. Additionally, the 
sputtering theory could explain the diffuse spectrum as the 
receipt of particles sputtered from a non differentially 
charged surface. Without the charging, the ions are not 
focussed into the higher channels, and can appear as a broad 
low energy distribution. The third phenomena, the shadow 
peak, can be described by tvo mechanisms. The first is a 
tr evonless detector, not scanning in angle, and observing 
опе trajectory, i.e., one energy. Alternatively, the 


detector could be scanning laterally across an emitted beam, 


ГЕ 


as opposed to vertically, and thus be picking up the same 
energy at each angle. Figure 308 illustrates the geometry. 

The problem with sputtering as a source, wasem a e 
had never seemed to be able to provide an adeguate flux of 
ions. However, the only previous mention of sputtering had 
been prior to the advent of the Sigmund theory, and before 
it was acknowledged that the magnetospheric regions 
contained large numbers of oxygen ions {Ref.1}. Preliminary 
calculations at this point indicated that the £lux 90g сай 
oxygen sputtering was sufficient to account for the observed 
phenomena. It then remained to model the situation and 
determine more accurately the flux at the detectors due to 
sputtering. 

It was now necessary to evaluate the physical 
characteristics of the environmental interaction with the 
Spacecraft. As the bulk of the vehicle surface was covered 
with solar panels, the most likely target of the incident 
ions would be amorphous silicon dioxide. For a magnetically 
disturbed time period, the ion flux would be composed 
primarily of oxygen and hydrogen. The Sigmund theory has a 
strong dependence on the mass of the incident ions, and it 
is reasonable to expect a much larger sputtering yield for 
the heavier ion. See, for example, Figure 14 for Н" апа 0" 
On amorphous silicon. 

The next step was to use the equations (1-7) and (1- 


8) to calculate the yield. One roadblock to this was the 
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Detector 
Scaming Vertically Through Energy 


Detector Scanning Across Angle Receiving 
Constant Energy 





Figure 38. Detector Beam Viewing Geometry 
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inability of the Sigmund theory to adequately treat 
sputtering from multi-atonic target: There are some data 
points available for impacting molecular oxygen от БИН 
and silicon dioxide, and they are listed in Тарғын 
оре. 3516 If the data from Table 7 are applied to Table 9, 
we see that the sputtering yield for O* on ві11соп ин 
range from .l to .18 at 10 keV. Similarly, the yield on 
silicon dioxide will range from .23 to .44 atoms perma 
10 keV. The sputtering yield calculated from the Sigmund 
theory for atomic oxygen on silicon is too high, oan x 
accurately approximates the yield from 510;. We will, 
however, take the worst case approach, and calculate the 
yield for silicon. To do this, we will increase the binding 
energy, to give us values more in line with those from the 
experiment. The two yield curves are displayed in Figure 
31. The new binding energy is 26.8 eV. 

Given that this new curve now reasonably describes 
the sputtering yield, the problem is to calculate how many 
particles emitted from a surface will reach a differentially 
charged detector for a realistic incident flux. The total 
yield of sputtered particles will be described by 
o) у = Ы SV эш а 27 

г 
In this expression, Г 15 the maximum energy a sputtered 
particle can have, that is, the maximum amount of energy 


transferred. This maximum is calculated by 
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Figure 31. 
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Г = 4 Mı Mz Em, 
(Мама) 

and is equivalent to the maximum energy transferred in а 
binary collision: The function £(v) is the distr 2 0 
function of the incident ions, and S(v,8) is the sp 2 
yield. 

With the total yield in hand, we may now attempt to 
calculate the countrate at the detector. The total number 
ОЕ counts the detector receives may be described by the 


following equation, 
(3 2) соната За" 


This assumes that the detector is differential in time, 
area, solid angle, and energy. If we integrate over the 


area increment and time, we get, 

(220) Counts=T дА шашы ыы, 

Thus, the’ countzarefis 

(3-4) Countrate=Counts/Tine=0\( (4 ЧАДА ШЫ 


To continue, we must assume that limits of integration on dR 
are not a function of energy, and that the value of the 
integral may be given by a simple number, characteristic of 
the detector. This is not actually true, as the limits on 
integration depend strongly on energy, since f(v) may not 


fill the entire 5 by 7 degree aperture of the detector. The 
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inaccuracy of this assumption vill surface later. With the 


assumption, we derive, 
E 5) Countrate=CR=6A бя gon 


The two constants, ôA and да, are combined into the 
geometric factor G, a constant for the detector. For 01171 
detectors, С=1.6 * 198^? m* ster. Due to the differential 
aperture of the detector, the value of the dot product ў- А 
вонеюргохітасе|у v (1.e. cos 921). 

We must now determine the distribution function of 
the emitted particles. From the earlier discussion of the 
ünmompson Theory, itis evident that the energy distribution 
Srecie yield flux is proportional Со E/(E+Ua)?. That is, we 


have the function 


$ ) ЧУ = А Е 


Е (E+Uəe ) 2 





ev 


me Elux Y, is also given by 
3-7) У з В le E £(E) dE 
e 


ЕГЕ F(E) is the distribution function as a function of 
energy of the emitted particles, and B is a constant that 
involves the integral over the solid angle. IF we combine 
these two relationships, it is evident (disregarding 


constants) that 


65-8) Е(Е) = ау - | 


1. 
Е 


Си 
СУ 
СУ 
+ 
С 
5 
W 


ү 


че will now determine the appropriate normalization 
constants. From the definition of eane distribution 
function; 


22208 n = ЈЕ Bey) о 


Continuing the calculation, we assume that the yield flux 
and рц distribution function have similar angular 
dependencies. Thus, £(v,9)=£(v)cos®. If we break dovme isa 


and substitute for energy, we get 


(3 10) ш: со58 сой. ОМ 2^ 
сок 

In the case of emission from a surface, © will range from 0 
Lo ins 3m O will vary from @ to т/2 Integrating over 
these limits, the equation reduces to 
(23:11) п-Ат(2/07)“ Ea dE 

(Bot Ua) 
This integral will extend the full energy range of the 
emitted particles, that is, from ð energy to r, the maximum 
energy transferred in the collisions. If one then 
substitutes a dummy variable, x*, for E the integral has a 
standard form (Бер 3/ colving, ме me m gele 


n=An(2/m? ) * Г jode: dE 
o (E + Ua)?” 


ВА (Мо ЇЕ хо ах s = г“ 
o (x7? + U)? 
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Ша ОООО Wwei(x- + Ја)7) t x/(8Use(x* * Ue)?) 


22750 289 шаш (х/ Џеј“) рэб 


эгэл UT) L/U = 2/(Џа + г)) 
НО а > A tan CrZUgj'*] 


Defining ß as the large constant on the right side, we get 


Ам)“ п 


OT 


25-12) nm <p 
e eng this normalization constant, the distribution function 


becomes 


(3-13) ЕСЫ, Ө) =жқп m?(2/n08X4)€cos0 1 
Rr Ua) 


If we then calculate the yield flux, Y, we obtain 


V (1/2m)% (n/B) IE E dE 
e (E t Ug)? 


№14) МЭ Үг 1221-01-51, /2172(6Г7 = Џеј“ + (17205) 3 


ПИО (n/B) 


Or 


8-15) пакету ВУ ó 


Nerefore, substituting in the distribution function, 


lo) AE Ө пе Y 1 со58 


2m 6 (E + Ue)? 
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We now have an expression for the distribution 


function of the emitted particles. This expression Gamma. 
be inserted with one modification into  egn. (3-5) 20115 
determine „спе countrate at the detectors. The necessary 


modification accounts for the differential charge between 
the emitting surface and the detectors. For the particles 
to be accelerated to the detector, the surface must be 
charged positively relative to the detectors, therefore Ф 
will be negative. Since E will map to E + qd (Lion ya 


Theorem), we will obtain 


22-21 ЕСЕР ЩЕ ПЕ со88 


2 6 ((E-%) + Ug)? 
Here, we have assumed that the ions are singly charged, and 
Ф will be converted to eV from У on this basis. The 


countrate expression is then, 


(3-18) CR=G(2/m?) | £(E,0,9) E dE 
(3-19) CR=2GY cos8 E dE 
го ((Е-Ф 12075 


To find the countrate, we must now integrate for each of 
the detector energy channels over its energy range. As this 
resolution is 205, the upper and lower limits are 1.1Е- апа 
0.9Е-, where Ee is the central energy of the energy channel. 
The integral in egn. (3-19) is a standard form, and the 


countrate may finally be written as, 
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(3-29) СЕЗ201/ 0059 21228 аа 9) + (1.ВЕ- + Џе - Ф) } 
тё ООО -—9j* -2029Е- + Uo - $)*? 


ШЕ спв energy channel limits are both less than 2, then the 
countrate over that channel is zero. ТЕ the lower limit is 
below 2, then the integral is taken between the upper limit 


апа Ф. 


Pee CALCULATION OF THE SPUTTERING FLUX AT THE DETECTOR 

Using the Sigmund Theory developed above, and the 
derived expression for the countrate at the detector, we 
developed a program to calculate the sputtering yield, 
incident and emitted flux, and the reponse of the detector 
to changes in both the energy channel and differential 
charge, 2. Data for the ambient plasma compositions was 
available from the Lockheed Ion Mass Spectrometer for days 
83, 86, апа 2909 of 1979. The values utilized are listed in 
Table 18. 

The first portion of the program vas the calculation of 
the total yield or emitted flux. This value is calculated 
From the following expression; 


(25-21) Ү з r оо ес”? у 


ә 
In this relation, v f(v) is the flux of the ambient ions as 
ction of velocity, апа S(v,8) is the sputtering yield 
as a function of velocity and angle. Thestactor Sn(e) іп 
S(v,8) exists in tabulated form, and is not described by а 


simple continuous function. If the points are fitted to a 
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curve, we arrive at the equation, 


(3=22) Snl€)= .35874 - (.17885)loge -(.14355° (Tos IR 
+(.06193)(loge)? + (.04228 Коле EP 
-(.00541) (loge) –(. 00630) 69а 
-(.00092)(1о46)” 

This eguation and the available tabulated points are 

displayed in Figure 32. Since S.(€) is a complicated 

function, ће integral was evaluated numerically. The 


function was calculated at a thousand points over an energy 
range that depended on the temperature of the ambient 
Maxwellian distribution. For a 4008 eV distribution, the 
high energy was 50000 eV and the low energy was 1080 eV. 100 
eV was chosen as the lower limit, as most sputtering yields 
are approximately zero by that point. At the 54900 eV 
роіпе, Епе value Of the distri burevon паре о која has reduced by 
a factor of 19-229 from the peak value. 

The yield obtained in this method was then multiplied by 


.02, as approximately 2%-4% of the emitted flux is  ionized 


ке Бони З ай This ion yield is substituted in the countrate 
integral to obtain the detector responses. The values of 
the integrals were varied over energy channel and 


differential charge to determine if the responses were very 
narrow as observed. Figures 33 shows the values of the 
incident, emitted, and ion fluxes for oxygen as a function 
Ot energy for dav 96 Of QS SE The figures for days 83 and 
200 are similar. From this it is very clear that the 


incident flux is much greater than the emitted ion | 
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Day 86 of 1979 


84 


Thus there is a net accumulation of positive charge on the 
Sipacecraft due to the ambient ion current. 
вашыя memo e тор ауо= о the countrate as a 

Fiction of energy channel for a constant шэ eV) 
differential charge, for days 86 and 200 for oxygen induced 
emission. We see that the countrate is high in 3-4 of the 
energy channels. The data indicates that it in fact should 
be concentrated in only one channel. However, the above 
equations assume that the detector surface is a flat plate, 
and that all particles reaching that plate are counted. 
foes is in fact not true. The detector viewing cone is 
probably not completely filled, as noted previously, and we 
are seeing the error due to treating the factor ба as а 
eamstant in eqn. (3-3). This factor will serve to reduce 
the countrate overall, and since particles at the edges of 
the viewing cone (with different energies from the cone 
center) will be reduced, the peak will most likely become 
narrower. The response for a specific channel and varying 
differential charges is quite similar. That is, a given 
channel will only respond for a certain narrow range of 
differential charge. Similar calculations for the ambient 
hydrogen population on the same days indicate that it can 
also generate a significant flux at the detectors, simply 
due EO the large amount present. The relative 
contributions of these two ambient plasma components very 


widely, depending on the temperature and density of their 
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Figure 35. 


Day 200, 1979 


O* Incident, 


87 


respective distribuel.ons, In general, the oxygen 
contribution is dominant. 

The countrates from the detectors are available for the 
above three days. On day 83 ot 1981, the mess 69 
countrates varied from 29 - 880 counts per second, with the 
average at 368. For day 86, the range was 24 - 1480, with 
the average at 3080. Day 248 was much lower, between 25 and 
213, with an average value of 124. The available mâss 
spectrometer data was averaged over one hour periods, 50 ме 
will generally calculate only one representative number for 
the charging time frame. The calculated countrates were 
380, 325, апа 400 for days 83, 86, and 290, respectively. 
Considering the approximate nature of the calculations, the 
theoretical responses compare favorably with the measured 
data, with the exception of day 290. However, it has 
already been noted that this data from the Hi detector 
indicated a sputtering source different from those seen on 
the Lo detector. This further supports that observation, as 
the sputtering yield for this event is clearly much lower. 

As the assumptions in this calculation tend to minimize 
the number of detected particles, it seems reasonable to 
conclude that sputtering from a differentially слац — 
surface is sufficient to account for the narrow triangle 
peaks. These discrete responses will serve to explain the 
triangle and shadow peaks. This leaves the diffuse 


background. We require a broad response, with the particles 
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scattered equally across the energy channels. Figure 36 
vr eserine <rogwdistrıbuclon function asa function of energy 
for ап emission angle of 0 degrees. тре stribution 
function has its maximum in the low energy channels, and 
numerically compares favorably with Figures 21 and 24. The 
Droad spectrum may then be described as the returning 
emissions from non-differentially charged surfaces. 

We must now explain the shadow peaks. These appear to be 
a combination of two occurrences, acceleration from a 
differentially charged surface to a motionless detector, or 
acceleration from a differentially charged surface to a 
detector scanning laterally across the emitted beam, as 
opposed to vertically. See Figure 31 for a view of the 
geometry. Since this is then a variation of the triangle 
peaks, the source and results would be the sane. 

The trajectories to the detectors can be recreated by 
the NASCAP program. Figure 37 displays the particle 
trajectories for the observed angles and energies at the Lo 
detector on day 86. From top to bottom, the lines are as 
follows; a)84 eV, 145 degrees, or 12.5 degrees above the 
spacecraft plane, Б) 70 еУ, 152 degrees, с) 60 eV, 175 
degrees, d) and 53 eV, 188 degrees. One can see that the 


particles track back to the spacecraft surface to different 


po nts. These points are all located in the зојак array 
portion of the spacecraft surface. These results are as 
expected. The solar array has charged differentially, and 
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Figure 37. NASCAP Trajectories, Day 86 of 1979 
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the sputtered particles folleowe specific. trajectories еса 
electric fields. The limited view cone of the десессо 
picks up specific energies. 

NASCAP was also run for the potentials observed on day 
94. The code was run to simulate the responses of the Lo, 
Hi, and Fix detectors to a range of incoming расата 
energies. The energies ranged from 2.0 to 90.1 8 Às 
noted above in the observations, the spacecraft vas charged 
to -60 V, and peaks vere noted at various lower energies on 
the detectors. The Hi detector vas parked at 67 degrees, 
and the Lo at 91 degrees. Thus, Hi and Lo detectors were 
sampling approximately the same region of space. 


Figures 38 and 39 depict the Lo detector trajectories in 


two different planes. The first shows the top of the 
spacecraft surface. Particles at energies less than 43 eV 
are returned to the spacecraft. The trajectories 


corresponding to observed particles are marked with an 
asterisk. The second figure is a complementing view from 
the Jett side of the first o figure. These two figures 


indicate that the actually observed particles are emitted 


from the satellite's forward end. Additionally, there are 
gaps from which no particles are seen at the detector. An 
explanation for this is the existence of patches of 


different nonconducting material on the forward surface for 
use in the Satellite Surface Potential Monitor (SSPM). When 


these patches (kapton, silvered teflon, and guartz) charge 
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Figure 38. NASCAP Trajectory, Day 94, Lo Detector 
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Figure 39. NASCAP Trajectory, Day 94, Lo Detector 
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differentially, tche particles reaching the detector from 
each will follow different trajectories, and have different 
energies. Additionally, the S-band antenna mast on the 
forward end is nonconducting, and may be a particle source. 

Figures 40 and 41 are similar views for a model of the 
Hi detector trajectories. We see that the origin of the 
particles is also on the upper surface, and the data shows 
the same striations. Again, the only nonconducting surfaces 
on the forward end are the SSPM and the antenna. The detail 
on NASCAP is not sufficient to study the extremely strong 
local electric fields expected in this region. If the 
satellite potential 15 190 У, then the field due to the 
sheath is approximately 1 V/m. In the forward surface 
vicinity, with a differential potential of 30 V, the fields 
will be on the order of 39-69 V/m, and will significantly 
Mmeseort the local trajectories of particles. Thus, the 
NASCAP trajectory plotting routine is somewhat unreliable in 
mars region. Looking at both figures, however, the origins 
for the two detectors are reasonably close, and allowing for 
the local fields, may possibly come from the same surface. 
As a working hypothesis, we will assume that the SSPM is the 
source of the emitted particles. 

Figures 42 and 43 model the fixed detector. The views 
are similar to those above. We see that the particles 
viewed by the fix detector come from completely different 


regions of the spacecraft. All the obviously  reattracted 
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Figure 40. МАЗСАР Trajectory, Day 94, Hi Detector 


96 


3.3 


31 


29 


211 


25 


23 


21 


197 


VE 


135i 


зш 


11 


єз 
“о 


л нч 94 07 





ПРО НЕ о у а ЗОРЕ] 12723 14 15 1617 
X-AXIS 


Figure 41. NASCAP Trajectory, Day 94, Hi Detector 
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Figure 42. NASCAP Trajectory, Day 94, Fix Detector 
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particles have energies of less than 43 eV, as did the other 
two detectorsk We can see that the trajectories о - 
particles skip over the bellyband region. This is 
reasonable, as the bellyband would be expected to be at the 
satellite mainframe potential, and thus particles emitted 
would not be accelerated by a differential charge. The 
trajectory calculations for the Fix detector also spp 
the theory ot sputtering trom different o nS charged 


surfaces. 
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IV. CONCLUSIONS 


It is reasonable to conclude that sputtering is a likely 
source of the observed phenomena. It can satisfactorily 


r plain the flux and Structure оғ the triangle and shadow 


peaks. It also can account for the observed broad spectrum 
e ons. These fluxes are large enough to be measured by 
the detectors, but are they large enough to cause 


substantial damage to the spacecraft surface? 

For tnis data set, the sputtering yield was less than 
the incident flux, so that the net result was particle 
accretion on the surface. Coincidentally, the data in this 
thesis were taken from relatively low energy days. On 
higher energy days, the sputtering yield increases, and the 
net result is surface depletion. However, the surface 
damage caused by this process appears to be minor. Assuming 
a higher energy day, with emission on the order of twice the 
incident flux, one atomic layer will be removed from a 
Silicon dioxide surface every 8 years. Since the solar 
cells are covered by the glass shields, they will be 
impervious to the surface sputtering. 

Another facet of this is the damage to surface 
conducting coatings on spacecraft. These coatings are only 
several atoms thick and may be expected to sputter from the 


Surface within several months of the satellite launch. This 
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may account for the observed charging оп = І 5 ни 
satellite. 

Since the energy Spectrum of the sputtered "parnm 
dies off quite rapidly, and the observed differential 
charges are usually not large, it is not likely chaca 
great contamination of the charging peak results. The 
sputtering in general will probably be uniform over a given 
material. Surfaces will be worn down over large areas as 
opposed to experiencing significant localized damage. 
Therefore, it appears that no short term damage vill result 
from sputtering on spacecraft. 

Further work is required to accurately calculate the 
sputtering yields from various surfaces to prove or Qi 
these conclusions. At this point, given the current state 
of the sputtering field, it is not possible to state with 
any certainty whether extensive damage will result from 
ambient ion bombardment. There appears to be no work in 
progress anywhere to determine experimentally the sputtering 
Yields for most ion-target combinations. Since the Sigmund 
Theory is so limited, these results are necessary for any 
actual calculations of the effect of the ambient heavy ion 
bombardment on space vehicles. 

This thesis demonstrates that sputtering is the probable 
cause of the observed low energy ion fluxes. Additional 
work is required to refine the NASCAP trajectory information 


and couple it with the theory developed here to model the 


1@2 


detector response accurately. Ber recommended that this 
(ве сопегпцесшипеп апа 1: Zurther applicable sputtering 


data become available. 
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APPENDIX 


TABLE 1 


MAGNETOSPHERIC IONS AND THEIR SOURCES {Ret oa 


Origin Source Tons 

sun Solar Wind H* Не Сы 

Earth Ionosphere Ht Het iog 

Earth High Altitude Н",Не",Не"" Ge] 
Thermal Plasma, 0575 ВЕ 


Plasmasphere 
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STOR MEME 


Date 


H* Flux 
ИЕ. > 


Density 
vd eem > 


Н” 
Не“ 
О” 


Mean Energy 
кеў 


H* 
Не" 


а 


RINON C URRENTECOMPOSITION iRef. 


a7 July 
1977 


12725 


11 


2072 


26 Моу 


ТАВЬЕ 2 


ПЕ 


4. 


44 


gu 


1795 


20 


09 Mar 


1978 


125 


143 


lU 


200 


0 


8) 


30 Мау 
12388 


6.0 


58 
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TABLE 3 


STORM TIME MAGNETOSPHERE COMPOSITION {Ref. 8) 


Date 04 May оны 29 Jan 82 June 
12977) 1977 1978 1978 

H* Flux e m 2.4 19 2:9 

Lo cemas — 

Density 

ГИ em З 
Hs 26 20 19 21! 
He 25 216 9125 231 
Нег” 228 > S 19 
Gi 340 ШЕ 4.9 60 

Mean Energy 

keV 
Бе gb 91.5 4.2 68 
мени PASE ЕН У 8.6 14.9 
Не“ 4.8 2. 6.4 4.8 
07 4.9 4.7 93,24 421 
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TABLE 4 


VALUES OF THE REDUCED STOPPING POWER {Ref.27} 


SU, © ЭБЕ) £ 
0.120 0.002 0.403 0.20 
0.154 0.004 0.405 0.40 
22211 0.01 0.356 1.00 
G2 on: 0.02 291 2.00 
даз 0.04 0.214 4.00 
2.372 0.10 2.128 10.0 
02395 D. us 0.0813 20.0 
0.0493 40.0 
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TABLE > 


CALCULATED AND MEASURED SPUTTERING YIELDS POR ЧАБТОВ 
MATERIALS {iRef. 23} 

THE NUMBERS IN PARENTHESES ARE THE YIELDS CALCULATED 
FROM THE STSHUNBZETREORT 


TARGET Sputtering Ratio (atoms сн. 
MATERIAL Ue(eV) Ze Ме“ Ar* KE Хе” 
Pb 2201 82 25 10-5 24.0 44.5 
222212 (44.0) (74.6) 

Ag 2.94 47 4.5 10.8 2305 36712 
20223 (20.1) (24.9) 

501 oan 50 5 4,3 Bub 11188 
(522250 (190257) (24722 

Cu 32246 29 9102 6.8 eo e 
СЭ 7) 22020 СІЗ ЗЕ (15250 

Ац Зао 19 5-6 20222 24:25 39 “И 
И) (23230 (2754 

ра 3 87 46 2223 53 1025 14.4 
(92223) (14.1) (1812 

Ке 4.29 26 Е x 4.0 4.9 
> 15) (32224 (1140 0 

Ni 4.43 28 1.4 E 5 7.6 
ОЕ Эл (9.01, (11728 

У DENS 23 0.3 1.@ I 105 
(3.285) 62 1) (8.9) 

РЕ 52132 78 1.9 S З 14153 16.90 
(9:39) (152488 (197128 

Мо 6.82 42 8.6 185 20 328 
5.24) (T DE @ 9; 

Ta 8.06 T3 8.7 16 3.1 4.8 
(6.62) (19.9) (12. 1 

Я 87710 74 1-0 253 4.7 6.4 
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TABLE © 


CALCULATED SPUTTERING YIELDS USING A 
MODIFIED BINDING ENERGY. 


TARGET Ре и Ratio (Atoms/ Ion 
MATERIAL 1.7 Џе Ат“ Kr Xe* 
Pb 3.42 1725 2959 43.9 
Ag 5240 7.47 ME 8 14.1 
Sn >: 29 ЕЕ 5 17285 24.1 
eu 5.88 294 7.84 СЕС ? 
Аџ 6.44 9.47 14.1 16.35 
Pd 6.58 5.44 9545 10.7 
Fe 45423 3-3 5.29 51:29 
Ni 919 SI 5229 6.94 
У 9.01 2x 25 3465 5.24 
РЕ oo 5-50 9.00 11.24 
Mo 1129 2242 4.53 5.94 
Ta 1.325 34:83 5.85 7.24 
W 14.79 2.99 5.48 7.99 
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TABLE 7 


RATIO OF SPUTTERING YIELD PER ATOM FOR MOLECULAR 
AND ATOMIC ION BOMBARDMENT (REF. 31] 


Projectiles S Ag Au 
CIT бй, 1.09 

Se-Sez 7225 1.44 1.44 
Te-Tez 1530 1 „82 2.02 
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TABLE 8 


SCATHA PARTICLE DETECTOR ENERGY CHANNELS {Ref. 4] 
Lo Detector 
Step Energy(eV) Step Energy(eV) 
@ -0.334 B 251225 
1 02229 33 30.10 
2 219592 3 34 34.57 
3 = ПЕС >> 39.10 
E -0.08 36 45.55 
5 0.00 37 ©2022 
6 0.09 38 60.00 
7 022 39 68.77 
8 0.34 40 18,55 
9 0.50 41 90.10 
10 @.67 42 10542552 
2-1 0.87 43 1192532 
12 22210 4 4 132.255 
148: 28:36 45 155.43 
14 12766 46 173.09 
15 2.00 47 203.98 
16 2.39 48 233292 
17 2.83 49 26720 
18 3. 34 50 306.30 
19 3.94 51 350.74 
20 4.60 52 401.96 
21 517 53 4579 3 
22 625 54 225 50 
23 725 55 601.05 
24 8.40 56 687571 
25 2) 2 57 7837570 
26 11823 58 992.13 
27 121595 59 1.039223 
28 14.96 60 1184.233 
29 17023 51 10255142 
30 19:83 62 15290585 
D 22.80 63 1777.60 


II 


Step 


von a CJ) NH ә 


TABLE 8 


CONTINUED 


Hi Detector 


Energy (eV) 


AR 
= 


El? 


Energy (eV) 
12052 
13855 
158659 
1820. 
2086. 
23328 
27427 
31:43: 
36017 
41278 
47287 
Sd 
6206. 
71888 
81428 
93252 
12681. 
12227 
14010. 
16045. 
183742 
21041. 
24096. 
27995» 
3159 WP 
361817 
414317 
47441. 
543238 
622034 
712265: 
315558 
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TABLE, 9 


THE SPUTTERING YIELDS OF 18 keV ION BEAMS ON 
ЗРЕТСОМ АМО 5 ПЕЕ ЗШ DISRIDE (Кет. 36) 


Target Incident Ion Beam 
Art O2* 
S 0.94 0.20 
SiO2 @.69 @.48 


ШЕКЕ 


TABLE 19 


AMBIENT PLASMA VALUES FOR DAYS 83, 86, AND 200 OF ||| и 
DATA FROM THE LOCKHEED ION MASS SPECTROMETER. 


Day Particle Number Density (m`?) Energy(keV) 
83/ H+ 3.845 * JP 8.852 
ШУ На” АВТ ЕКИ 129.778 
Не 2.865 Lg 33 ЗЕН 
O 3.188 * ие 4.21 
86/ Н 2.606 X 77023 5.885 
12973 На“ 2.850 *' 105 43.728 
He 2.418 * 102 3.662 
Оз 1,635 x (IOS 4.035 
200/ Н+ 9.064 * 109 5.278 
1979 Ha* 4.209. * 105 1.273 
He 2.1697 % 7192 0.420 
бум 4.093 * 10° 1.957 
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